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It is shown that a Bose gas of two-level atoms in the intense resonant laser field at zero 
temperature is a mixture of two condensates with a definite ratio of the densities. The criteria 
of stability are found for the stationary states of such system against the increment of the 
amplitudes of quasi Bogoliubov elementary excitations. Besides the usual acoustic mode the 
gap mode is shown to exist and the magnitude of the gap is proportional to the laser field 
amplitude. The involvement of the gas nonideality under definite conditions results in an 
instability and decay of the condensates. 

PACS: 03.75.Kk, 03.75.Mn 
1. Introduction 

After the experimental realization of the Bose-Einstein condensation in a gas Q| and later in 
a binary mixture of Bose gases P| in magnetic traps at ultralow temperatures the adequate the- 
oretical description of such systems is undoubtedly actual. This is confirmed by a great number 
of theoretical papers on the classification lj| of all states of a Bose mixture (within the Thomas- 
Fermi approximation) and its dynamics (longwave collective excitations |J|, metastable states 
the spatial separation of the mixture components In most papers the mixture is described 

by two Gross-Pitaevskii equations and the concrete distinction between the mixture components is 
not considered. The representation of the mixture components as atoms in theground and excited 
states is employed in the works on studying the scattering ^ and absorption {J of the laser radia- 
tion by the Bose condensate. However, the question about an existence of the stationary coherent 
state of the interacting system, "mixture -|- laser field", has not been considered yet. The present 
work is devoted to the theoretical investigation of this problem. In particular, the equilibrium ratio 
of the condensate densities is determined as a function of the resonant laser field amplitude and 
proximity to the saturation, the spectrum of quasi Bogoliubov elementary excitations is found, and 
the criteria of stability in the system of two condensates in the laser field are determined with the 
involvement of the recoil momentum and nonideality of a gas. 
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Let us consider a weakly non-ideal Bose gas of two-level atoms in the resonant laser field of the 
high intensity at zero temperature. The monochromatic laser field, acting on the Bose condensate 
of atoms in the ground state (for simplicity, we neglect noncondensate particles), converts some 
fraction of the atoms into the dipole-excited state so that each excited atom moves with the same 
velocity governed by the recoil momentum. If the laser field intensity is so high that one can 
neglect the spontaneous decay of excited atoms with respect to the induced emission, then, during 
the small time as compared with the time between collisions of excited atoms with quasiparticles 
of the condensate of unexcited atoms, the macroscopically large number of dipole-excited atoms 
moving with the same velocity is formed and the second Bose condensate appears. Below we 
do not treat the process of the formation of this condensate and assume that the state with two 
condensate has been already formed, taking as a condition that the natural linewidth of dipole 
transition is 7 0. The criterion of a smallness for 7 will be given in Sec. 2. At last, from the 
viewpoint of maximal simplicity of the effect treatment, the atomic system is considered as spatially 
homogeneous. 

2. The ground state of two-level gas in the laser field 

The Hamiltonian for the system consisting of Na unexcited atoms and A^^ dipole-excited atoms 
in the laser field {V = l,h = 1) reads 

Htot = Hgas + Hint + Hpfi, Hgas = Ha + Hf, + Hab, 

P Pl+P2=P3+P4 
P Pl+P2=P3+P4 

Hab = E Uaba^^b+^bp^Qp^. 

Pl+P2=P3+P4 

Here Op and bp are the atom annihilation operators with momentum p in the ground and dipole- 
excited states, respectively, and Uik = Airaik/m are the parameters of the pair interaction expressed 
in terms of the corresponding scattering lengths. To exclude the problem of collapse, we suppose 
J7jfe > 0. The interaction of the gas with a single-mode laser field is described by the Hamiltonian 

Hint = ^{gkcla^bp+k + glbp^k^pCk), (2) 
p 

where and are the absorption and production operators of a photon with wave vector k and 

frequency ujk- For simplicity, the laser radiation is assumed to be Unearly polarized in the direction 

of unit polarization vector e^. The atom-field coupling constant has a usual form 

2 ^ ^ 

\gk\ = ^-KUk {ek,d) 
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and is related to the natural linewidth of the excited level for a single atom via dipole matrix 
element d, 7 = ^kQcP where ko = To simplify the notation, the factor exp(iA:r) in the matrix 
element of the interaction between a single atom and laser field is omitted. Lastly, the Hamiltonian 
of free photons is given by 

Hph = LJkC^Ck- 

The dynamics of the system is described by the equations 

^"g^ = ~^ X] UaaO-p^Clp^ap^ + ^ Uabbp-^bp^^ttp^ -\- QkC^bp+ki 

V+P2=VS+PA Pi+P=P3+Pi 

^^ = (^+^o)&p+ E Ubbbj^bp^bp^ + E Uabaj-,bp^ap^ + glap^kCk, (3) 

P+P2=P3+Pi Pl+P=P3+P4 

= ^kCk + J2 9ka'^bp+k- 
p 

Let us consider the system which consists of three coherent subsystems, namely, two condensates 
produced by atoms in the ground and dipole-excited states and the laser field. Let us find the 
ground state of this system, neglecting noncondensate particles. In what follows, it is more conve- 
nient to deal with usual complex quantities rather than with operators. Therefore the use of the 
representation of coherent states is appropriate. Let us determine a coherent state as 

\a) = exp(^apa+) |0) , 
p 

which can be considered as a direct product of all eigenvectors \ap) so that ap\ap) = 
Up \ap) , (opl dp = (opl a* and concerning the states with different momenta. Vector (|0) denotes a 
direct product of all empty atom and photon states. For dipole-excited atoms and photons, the 
coherent states are determined similarly: 

|6) = eMYl ^pK) |0) ' 1^=) = exp(cfec+) |0) . 
p 

Taking the recoil into account, condensates of atoms in the ground and dipole-excited states should 
move relatively each other. Let us denote the momenta of atoms belonging to these condensates pa 
and pb = Pa+k, respectively. If the gas is at rest as a whole, pa = —kNi,/{Na+Nb). Postmultiplying 
each of equations ((SI by vector \a) \b) |c), we obtain the set of three equations for the complex fields 

^Pa 1 bpi, , Cfc : 

= (fi + UaaNa + UabN,)ap^ + QkClbp,, 

= + ^0 + UbbNb + UabNa)bp, + glap^Ck, (4) 
= ^kCk +9ka*pJ)pt,. 
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Hereafter we consider a gas in the condensate state. Therefore Na and Nh are macroscopically large 
quantities, i.e., Na, Nf, » 1. We seek a stationary solution in the form 



an„ = ae 



bp. = be-''^\ Ck = ce-'^K (5) 



Pa — ) "Pb 

Note that a stationary solution (O takes place only under condition of the exact resonance between 
the renormalized laser frequency and the difference in energies of the excited and ground states of 
an atom 

W = - Ea- (6) 

From we get 

Eaa = + gkC*b, 

Ebb = Ubb + glac, (7) 
{Eb - Ea)c = LOkC + gka*b. 

Here the following notations are introduced 

^a = ^ + UaaNa + U^bNb, Qb = ^ + UJo + UabNa + UbbNb, Na = \af , Nb = \bf . 

2m Am 

Quantities Ea, Eb and Na (or Nb since we consider that the given total gas density is N = Na + Nb = 
const) are unknown. The third equation in ^ determines the renormalization of the laser frequency 
due to dipole transitions. This renormalization can be neglected if the laser field intensity is 
sufficiently high (formally, this means that Hp^ ^ Hint, implying matrix elements) so that 

1\T I |2 ^ \9k\'^ NaNb 

Nc = \c\ > 2 ■ (^) 

The condition of the smallness of the natural linewidth results directly from this inequality 

Nek?, 

Inequality ^ is worthwhile to compare with the initially implied condition of the sufficiently high 
field intensity in order to neglect spontaneous decays with respect to the induced transitions. In 
the terms of the notations introduced this condition is given by 

Nc > -k^. (10) 
vr 

Note that, if one takes resonance approximation {\ujk — ujq\ ^ uJk,^^o) into account, inequality (jTIHl 
yields the condition known in the laser physics for saturation of the population of excited state, 
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when an atom with the equal probabihty is in the ground and excited states so that Na = = N/2. 
The inequahty is expressed quantitatively as 



\uJk - wol ~ 7 < 1^1 



(11) 



where \A\ = \gk\ is the Rabi frequency in the conventional terminology. From (jTTIl we obtain 



(12) 



From the comparison of conditions Q and ifT^ one can see that they are equivalent to each other 
if kQ > i.e., when the volume per one atom is larger than Aq. According to (dU, this means 
qualitatively that a single atom is fallen at many photons. Note, however, that condition (fTTjl looses 
its meaning if the renormalization of the transition frequency is larger than the natural linewidth 
for a single atom, |ef, — — ujq\ > 7. In this case, it seems that 7 should be regarded as a total 
natural linewidth with respect to the collisional broadening. The calculation of the corresponding 
contribution is made in j^. 

In what follows, we assume that inequality ^ is fulfilled and Uk- Denoting e = Sa and 
re-denoting fife — > 0^ — 0;, we have from the equations (0) 

{e-na){e-VL,) = \A\' . 



Hence 



1 

'=2 



,7? = ±1. 



(13) 



Here the quantity 



(5a; = — = — Wo + 



e fNh-Na 



2m \Na + Nb^ 

plays a role of the deviation from the full saturation. The relation between the amplitudes of the 
condensates is determined from Q 



5u) 
'2\A\ 



5uj 
2\A\ 



+ 1 



a. 



It follows from above that the system in the laser field is in one of the states with the fixed ratio 
of the condensate densities (at the given sign of product rjSuj) 

2 



6uj 

-V — TTT + 

2 A 



5uj 
2\A\ 



+ 1 



(14) 
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When Su! = 0, the resonance field frequency equals u> = ujq + ^{Ubb — Uaa)N. Thus, if Uaa — Ubb, it 
is natural to take difference to — uq as an external parameter characterizing the deviation from the 
saturation. Let us introduce the dimensionless variables 

and write equation ifTHl in the form 



l + y= (^-V{x + zy) + ^j{x + zyf + 1^ (1 - y). (16) 

Equation (|TH|l is a quartic one with respect to y{x). The solution relative to the inverse function 
X (y) has a simple form 

x = -(z + -J=)y. (17) 

yi - V 

It follows from the above expression that there exists a symmetry with respect to the simultaneous 
substitution x — > —x and y — > —y. This results from the symmetry of the Hamiltonian of the 
system with regard to substitution a ^ b. In its turn, such symmetry between "up" and "down" is 
a result of the assumption about the smallness of 7. Note that there is no symmetry in the solutions 
with regard to the variation of the sign r/. For the branch rj = +1, quantity y (x) is a one-to-one 
function of x and, for ry = — 1, function y (x) at z > 1 within the interval 

takes four values (see Fig.l). Accordingly, for the branch 77 = — 1 at z > 1 there exists limiting 



values ^(ibxiim) = =F\/ 1 — z with an anomalous S-shape of the line in comparison with the case 
z<l. 

For x = 0, both condensates have either the same density, y = 0, for any sign or 



y = ±Vl-z-2, (18) 

at r/ = —1. Thus, the satellite states are possible if the coupling between atoms and field is 
sufficiently small compared with the recoil energy Er = (z > 1, being z ^ gjAj ^"-^ ~ Uaa, Ubb)- 
If the coupling between atoms and field is strong, z < 1, there are solutions only in the vicinity of 
the centre (y = — j^x) at |x| <C 1 (i.e., in the case \eb — £a — "^ol < 7 |^|)- 

So, it is seen from the analysis of function x (y) that (i) under the condition of saturation 
|x| <C 1 there exist four stationary states, where in addition to the usual solutions y ^ there are 
two satellite states in which populations of the ground and excited states are essentially different. 
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(ii) Beyond the region of saturation at |x| > xum there are two stationary states, and at |x| < xum 
there are four ones again. It is also interesting to trace the behavior y{z) at re = 0. If z > 1, 
bifurcation, associated with a possibihty of satelhte states, appears in the plot y{z) in addition to 
y = (see Fig. 2). 

3. Combined oscillations of the condensates in the laser field 

The stationary state obtained is unstable if the spectrum of elementary excitations of the system 
in this state has an imaginary part. To clarify a role of the interaction between atoms, we start 
from the case of an ideal gas. In addition, we will consider the question about the stability in the 
case of full saturation u = ujq, Na = Nb = N/2 and in the case of a satellite state. 

First, we find a set of equations determining spectrum of single-particle excitations of the con- 
densates. Let us write equations (0 for the complex fields in the linear approximation in the 
amplitudes of noncondensate particles in the following representation (hereafter g is a quasiparticle 
momentum): 

ap = apexp{-iet) , p = pa±q, 

bp = bp exp {—iet — iujt) , p = pi,± q, 

in which the system (jSJ takes a form of a closed system with constant coefficients (vector notations 
are omitted): 



dapa+q _ I (Pa + g) 



2 



-^T- - ( 2m - ^ + ^Maa + Vah^h)ap,+q + ^aa^^-q + (^afe + ^) V+g + ^^-ahbp^.q, 



dbp^+q ^ / {Pb+qf 
dt y 2m 

^~K,-q _ f (Pb-q? 



-I 



e- {uJ-UJo)+ 2flbb + UabNa)bp^+q + fJ-bbbp^-q + (fJ-ab + ^)ap,+g + HabO-p^-q, 

'dr~ - \~^2Wi e-{u: -ujq) + 2nbb + UabNa)bpf^-q + fJ-bbbpb+q + f^abo-p^+q + {fJ-ab + ^)aX-g' 

(19) 



where the notations fiik = Uik^J NiNk and A = glc are used. The phase of laser field is chosen so 
that quantity A is real and positive. 

In the case of an ideal gas, when all Uik = 0, the equation for eigenvalues of the system ifT!?]! . 
with involving ifT^ and ifT^ . yields 



m / 
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Here 5uj = tu — ujq + -^y. For the full saturation {6uj = 0, y = 0), we have 



E 



1,2 



2m " ' ' 



+ \A\', S3^4 = -|-+ry|A|± 
2m 



( M. 



+ \A\ 



Hence, for the given rj there is a gap 2 \A\ for two modes in the excitation spectrum at = 0. The 
gapless modes at ^ = and q — > yield a sound dispersion law 



E{q) ~ CsQ, 



\kq\ . 



where Cg = is the sound velocity. Far from the saturation (|x| ^> 1, |y| ~ 1) in lfT5|l for the 
given rj the two branches have a gap of about |a; — wqI and the others are gapless as before. 

Thus, in the approximation of an ideal gas the stationary condensate state is always stable in 
the sense mentioned above. 

In the case of nonideal gas one can find an analytical expression for the spectrum of the system 
ifT!?]) under full saturation, assuming that Uaa = Uht- 



E± = (P' - + A{A + 2^,fe) + Q'±2yJ (p2 _ ^2)q2 + (^p + ^^^(p _ 

P=-^-r]A + fl, Q=^, jJi=\UaaN, jiah = \UabN. 

2m 2m 2 2 

2 

First of all, we consider the long wave limit q — > and ^ S^(0) + |;^A^ where 

El{Q) =e-f^'^ + A{A + 2fiab) ± le - WabI , ^ = fx-Ari, 
AE = 2^ + er cos2 e±( ^^^e^ cos^ + 2A{A + fiab) sign - rjfiab) ) • 



(20) 



(21) 



Here 9 is the angle between vectors k and q . The critical value, in the sense of stability, is E'^. In 
the case when the argument of modulus in ipTjl is positive fJ- > r] {fiab + A) we obtain 



E^^{0) = 2A{A- fi + fiab){l+v), 

AE = 2{fi- flab- A{1 + ??)) + er cos^ 



1 



(A - 2??/x) 



fJ'-V ifJ-ab + A)_ 

It is seen that there is an instability in the state = 1 at /i > fiab + A and in the state rj = —1 

AE = 2{fi-fiab){l-i:- 



IJ- + tJ-ab + A J ' 

at /i < flab- This coincides with the known condition of instability of a binary Bose-gas mixture 
against spatial separation of its components. The instability at Er cos^ 9 > 2{fi + jiab + A) means 
that for the strong recoil the instability develops at first in the direction of the recoil momentum. 
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The corresponding plot of the spectrum IpT]) at cos^ 9 = I is given in Fig. 3. The opposite case 
/i < 7/ {flab + A) has a sense only for = 1 and does not result in an instability 

E^{0)=0, Aj; = 2(^ + Mfl + ^ /:'°''^ ) >0. 

V 2A + l^ab-fJ-J 

Let us now consider the case of transverse excitations with respect to vector k when Q = 

2m 



^2 



For Et at {A + [lab) - V^) + Afi > 0, we get 
It is seen that in the state ry = — 1 an instability appears for the condition obtained above /x < fiab- 

2 

In the state rj = 1 there is a region of instability 2A — 2 (/i — fj,ab) < ^ < 2^ with the centre at 

2 

the point = 2A — [jj. — fiab) where i?^(Q'min) = — (/^ — /^ab) • In addition, for g = in this state 



2m 

?2 . 



there appears an instability \i ^> A + jjab (however, it is always i?i(gmin) < E_{0)). The typical 
plot of the spectrum for this case (PT)]! is given in Fig. 4. 



In the case {A + fiab) ( ^ — V^) + < we have 



2 

i^^ = (|^+A(l-.)) (|^+^(l-,)+2(^ + ;.., 

whence it follows that the both states rj = ±1 are stable. 

Finally, let us consider the stability of satellites at Uaa = t^tb- Far from the saturation 
{\lo — ujo\ ^ A, fi, flab) one can neglect the nonideality of the gas and, as we convinced above, 
the condensates are stable. In the region of saturation |x| <C 1, the satellites exist only in the 
state with r] = —1. For simplicity, we restrict ourselves with the case when the density of one of 
condensates is much less than that of the other. For definiteness, Nh <^ Na (y ~ — 1). It turns 
out that one can derive an analytical expression for the spectrum of elementary excitations in the 
approximation of heavy atoms and strong atom-atom coupling compared with the recoil, putting 
formally k <^ q. Note that in these conditions it is not allowed to employ the long wave limit 



?2 



El = ^ (w'^ + V^ + 2A^ -fi^± ^{W^ - _ ^2)2 ^ (^i^w + Vf - ) ) , (23) 




W = ^ + ^ (3// - Mab) + ^Jlitl-flabf+A^, 



^ = ^ - ^ (/^ - f^ab) + J\if^- l^abf + A^, 
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where ^ = Uaa^, tJ-ab = UabN. Rewriting expression in the form 



") 



one can straightforwardly indicate the region of possible negative values , 



Hence 



liV > {WV - A^) . 



(24) 



Note that the condition never holds true at ~ ^ab and the satellite state is stable. 
4. Conclusion 

The main result of the work is the determination of a stationary solution of the equations 
of motion for the united coherent system of two-level atom Bose condensates and laser field of 
high intensity. In addition, the dependence of the condensate densities on the Rabi frequency and 
deviation from the full saturation are found. The analysis is performed for the stability of the system 
against appearance of the imaginary part in the spectrum of elementary excitations. This results in 
an exponential increase of oscillations, in the heating of the system and decay of the condensates, 
or in the transition of the system into one of the stable states. The question on evolution of the 
system after the decay of an unstable stationary state is a subject of separate study. 

The condition of applicability of the theory in temperature is given by 



Note that in the satellite state, when the condensate densities differ significantly, the condensate of 
the lower density at T > Tc converts into noncondensate particles which are additional with respect 
to the condensate of higher density. The question whether this condensate will be destroyed with 
adding extrinsic noncondensate particles represents an independent interest. 

The authors are grateful to Yu. Kagan and S. N. Burmistrov for valuable comments. The work 
is supported by the grants of Russian Foundation of Basic Research and INT AS. 
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n = min {Na, N^) . 



m 



[1] M. H. Anderson, J. R. Ensher, M. R. Matthews et al, Science 269, 198 (1995); K. B. Davis, M. -O. 
Mewes, M. R. Andrews et al, Phys. Rev. Lett. 75, 3969 (1995); C. C. Bradley, C. A. Sackett, and R. 



11 



G. Hulet, Phys. Rev. Lett. 78, 985 (1997); D. G. Fried, Th. C. Killian, L. Willmann et al, Phys. Rev. 
Lett. 81, 3811 (1998). 

[2] D. S. Hall, M. R. Matthews, J. R. Ensher et al, Phys. Rev. Lett. 81, 1539 (1998). 

[3] Tin-Lin Ho and V. B. Shenoy, Phys. Rev. Lett. 77, 3276 (1996); M. Trippenbach, K. Goral, K. Rzazewski 

et al, J. Phys. B 33, 4017 (2000). 
[4] Th. Busch, J. L Cirac, V. M. Perez-Garcia and P. Zoller, Phys. Rev. A 56, 2978 (1997); R. Graham 

and D. Walls, Phys. Rev. A 57, 484 (1998). 
[5] H. Pu and N. P. Bigelow, Phys. Rev. Lett. 80, 1134 (1998). 

[6] P. Ao and S. T. Chui, Phys. Rev. A 58, 4836 (1998); E. Timmermans, Phys. Rev. Lett. 81, 5718 (1998); 

H. Shi, Wei-Mou Zheng and Siu-Tat Chui, Phys. Rev. A 61, 063613 (1998); R. A. Barankov, Phys. Rev. 
A 66, 013612 (2002). 

[7] J. Ruostekoski, M. J. Collett, R. Graham and D. Walls, Phys. Rev. A 57, 511 (1998). 
[8] M. O. Oktel, L. S. Levitov, Phys. Rev. Lett. 83, 6 (1999). 



12 



Figure Captions 

Fig.l. The dependence of the relative difference in populations, y, versus deviation x from the 
full saturation. 

Fig. 2. The dependence of the relative difference in populations, y, versus parameter z (see ifT^ V 

Fig. 3. The example of unstable spectrum in the long wave limit for the state r/ = — 1 under full 
saturation at e,. > 2 (/i + ^ab + A) = 0.436^, and /x : ^ab A = 1 : 0.5 : 20. The scale on the vertical 
axis for the curve ry = 1 is reduced by a factor of 40. 

Fig. 4. An example of the unstable spectrum for the transverse excitations in the state r/ = 1 
under full saturaion for Q = 0, A = 2.5£r, ■ : ^ = 3 : 1 : 5. The scale of the vertical axis of 
curve r/ = — 1 is reduced by a factor of 4. The cusp in the curve E'^, rj = I, corresponds to zero 
modulus in 
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